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ABSTRACT 



The Airborne Digital Data Acquisition System (ADDAS) was developed by the 
joint Naval Postgraduate School/Harvard Unversity research project, called the Ocean 
Prediction Through Observation, Modeling, and Analysis (OPTOMA) Program, for the 
real-timx synoptic mapping of oceanic features. 

Given the capabilities of such a system, an aerial sampling strategy was 
developed based on the spatial scales of oceanic vanables and the aircraft and .A.DDAS 
limitations. The strategy was applied to "simulated" AXBT survey flights on 
"synthetic" fields. The simulated AXBT observations were objectively analyzed and 
compared to the synthetic fieid. By varying the simulated AXBT spacing and the 
objective analysis parameters, an "optimum" AXBT spacing to adequately represent the 
synthetic structure of the domain was deduced. 

To further develop such a system for both airborne and shipboard use. a 
proposed Ship/ Aircraft Data .Acquisition. Display, and Analysis System (SAD.AD.AS) 
’.vas presented. The following synoptic mappings produced could provide a possible 
tactical decision aid for antisubmarine warfare (ASW) operations: 1) Sea surface 

temperature; 2) Temperature/Sound speed at selected depths; 3) Maps of mixed layer 
depth; 4) Subsurface ducts/sound channels; and 5) Vertical sections of 
temperature/ sound speed. 
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I. INTRODUCTION 



The development of an Ocean Descriplive/Predictive System (ODPS) for the 
study and forecasting of oceanic synoptic mesoscale features and the study of the 
California Current System are the goals of the Ocean Prediction Through Observation, 
Modeiing, and Analysis (OPTOMA) Program, a cooperative effort between the Naval 
Postgraduate School (NPS) and Harvard University. The Airborne Digital Data 
Acquisition System (ADDAS) was developed at NPS in response to the need for a 
real-time, large-scale synoptic mapping of the mesoscale features for both the 
initialization of dynamic model forecasts and the verification of the dynamic models 
used. 

A. OBJECTIVES 

To develop a system such as the .ADDAS for use in operational antisubmarine 
warfare fASW) is the primary objective of this thesis. With the capabilities of such a 
system: 1) What optimal aenai sampling strategy should be utilized to acquire oceanic 
data for naval planning and operations? Z) What funner enhancements could be 
developed for a system like the ADDAS, and how could such a system be applied to 
shipboard operations as well? 3) What operationally applicable synoptic fields can be 
developed from the data, and how could such fields be applied in planning and 
execution of Naval operations, specifically in ASW? 

B. BACKGROUND 

With the availability of expendable aircraft deployable instruments such as 
airborne expendable bathythermographs (AXBTs), sound speed profilers (AXSVPs), 
current profilers (AXCPs), and soon-to-be-available air-launched conductivity, 
temperature, depth profilers (AXCTDs), the aircraft is the platform of choice for near 
real-time data acquisition. The P-3 Orion ASW patrol aircraft supports a system such 
as the ADDAS exceptionally well. It has the payload capacity and avionics for air 
launched sensors as well as substantial on-station endurance and all-weather capability. 
It can cover a sizable oceanic region in a matter of hours as compared to days for a 
research vessel. The ADDAS has been deployed on six OPTOMA P-3 missions to date 
and is a proven (research) operational system. It has been invaluable as a primary 
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means of data acquisition for the OPTOMA Program. Additionally, carrier-based 
ASW aircraft such as the S-3 Viking, and rotary-wing ASW aircraft including the SH-2 
Sea Sprite, SH-3 Sea King, and the SH-60 Seahawk, can provide additional airborne 
sensor deployment platforms. 

C. MOTIVATION 

The capability for near real-time mapping of oceanic variables has a direct 
application to ASW operations. As submarines become quieter, the warfare 
commander will have an even more difficult task of localizing and tracking them. An 
increased knowledge of the submariners' environment uill be required. An AXBT 
survey flight can provide a "nowcast" of the region for use in tactical briefings, and 
such planning should be a key factor in ASW operations. With maps from such 
AXBT surveys, the crfects of the ocean structure on standard tactics can be addressed 
and appropriate modifications to tactical doctrine can be suggested. Additionally, 
subsequent survey flig.hts could be completed to identify the temporal and spatial 
changes over the course of the exercise or operation and to refresh the original survey. 

Traditionally. ASW tactics and training are based on a "generic" ocean structure, 
and usually only one or two AXBT's are deployed in an exercise which tends to 
reenforce this concept. One may argue that a single AXBT could be representative of 
the ocean thermal structure in some region. However, if an AXBT was deployed in the 
middle of an eddy or in an area of strong temperature gradients, as are observed in the 
Gulf Stream, it may not be representative of that region. In the research field, the 
ADDAS has proved itself as a powerful tool. It is time for this capability to be 
developed for operational use in the fleet. 

The current version of the ADDAS and similar systems in use today are 
presented, followed by the oceanographic application of the ADDAS in the OPTOMA 
Program. Next, an optimal aerial sampling strategy is developed and applied to 
"synthetic" oceanic domains and compared through use of the objective analysis 
technique. A proposed Ship/Aircraft Data Acquisition, Display, and Analysis System 
(SADADAS) concludes, with application to both land-based and shipboard aircraft 
operations. 
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II. THE AIRBORNE DIGITAL DATA ACQUISITION SYSTEM (ADDAS) 



A. ADDAS DEVELOPMENT 

The OPTOMA Program's Airborne Thermal Structure Mapping Project at NPS 
culminated •with the development of the ADDAS and its successful research use. In 
the early survey flights of this project, AXBT analog traces were recorded onboard the 
P-3 aircraft by the AQA-7 lofargram svstem recorder. The traces were hand digitized 
using gridded tempiates based on the standard Na’vp/ frequency-to-temperature ana 
fail-rate-to-depth equations. This produced traces of coarse resolution and resulted in 
a substantial time-lag for the generation of "real-time" synoptic maps. 

At the Scnpps Institute of Oceanograpny, Vlr. Mereaith H. Sessions aad 
developed a system for the recoraing ana real-time digital conversion or' an AX3T 
analog signal, ivhich had been used successfully in the North Pacific Experiment 
(NORPAX) (Sessions, 1980). The Scripps system was made available to the OPTOMA 
Program and was 'used for the post-ilight digitization of analog tape recordings of the 

A. X3T signals acquired by the OPTOMAliP fiight on 18 July 1984. Inspired by the 
relative ease of data acquisition and the improved resolution of AXBT profiles 
provided by the Scripps system, as well as the recent development of an AXBT 
digitizing circuit for the Sippican MK9 digitizing unit, a prototype ADDAS became 
feasible and was developed at NPS (Colton and Mooers, 1985). First deployed on the 
OPTOMA13P flight on 27 October 1984, the ADDAS has proven invaluable in 
acquiring real-time data for synoptic mappings and model initialization. 

B. THE PRESENT ADDAS 

The ADDAS is a self-contained, semi-portable data acquisition, processing, and 
analysis system. It utilizes off-the-shelf hardware and modified shipboard XBT and 
CTD data acquisition and analysis software developed by the OPTOMA Program at 
NPS. The heart of the system is a HP9816 MC68000 16-bit microprocessor with a 
HP9121 dual 3.5 inch microfloppy disk drive. Aircraft 28 VDC power is inverted to 
115 VAC, 60 Hz power by an internal static inverter. AXBT drop positions can be 
entered either manually by the operator through the computer keyboard (latitude and 
longitude, as determined by the Tactical Coordinator/Navigator from aircraft inertial 
navigation) or automatically by the ADDAS's optional LORAN C unit. The AXBT 
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audio signal is received through the ARR-72 aircraft sonobuoy receiver, digitized in 
real-time by a Sippican MK9 front-end processor, and the trace is displayed on the 
computer monitor. Once the AXBT has completed its fall, the digitized profile is 
stored on a data diskette. Currently, the ADDAS is capable of digitizing only one 
AXBT at a time. To simulate a multi-channel capability, one channel of a stereo hi-fi 
video cassette recorder (VCR) is used to record the simultaneous audio signal from a 
second AXBT (of different VHP carrier channel) for post-flight digitization and 
analysis. The audio signal of the digitized AXBT is also recorded on the other VCR 
channel as a backup. .An oscilloscope provides a visual display of the .AXBT audio 
signal to help detect .fiilures both in the aircraft sonobuoy receiving system and in the 
AXBTs. Once the survey is complete, maps of objectively analyzed variables are 
produced on the monitor and hard-copy output is generated by a HP7‘^70A two-pen 
plotter, either on the return flight or on the ground. .A list of .ADDAS hardware 
comoonents and dteir current retail cost is inciudea in Appendix A. 



C. ADDAS LIMITATIONS 

Limitations to the system can be classed into two categories, intrinsic and 
extrinsic. These limitations affect .ADDAS operation and ultimately the density of 
AXBT coverage in a surv'ey. 

1. Intrinsic Limitations 

These limitations are specifically inherent in the system design and operation. 
They are; 

• Not multi-channel capable (presently). Two AXBT's of different channels cannot 
be digitized simultaneously. An additional MK9 or the development of a 
multi-channel MK9 would be required. In addition, an advanced microprocessor 
with a multi-tasking capability is required. These are very feasible with today's 
technology. 

• No hardware redundancy. Although the system has proven exceptionally 
reliable, if one component should fail the real-time capability will be lost. 

• Inverter operation. The inverter is rated at 5.6 amps. With only the 
microprocessor, disk drive^ MK9, and VCR in the circuit, the inverter operates 
at 6.3 amps continuous. This has not been a problem to date since the .VCR is 
only used during data acquisition and it is disconnected when the plotter is being 
used. An inverter with a higher current rating may be required. 

2. Extrinsic Limitations 

These limitations affect the overall capability of the system to perform its 
given mission and are related to the sampling platform and the sensors employed. 
They are: 

® Aircraft weight and balance. Aircraft are constrained by a maximum gross 
weight for take off. It is a function of the environmental conditions at the 
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departure point. For navy aircraft, these limitations are .specified in the 
respective NATOPS (Naval Air Training and Operations Procedures 
Standardization) Manual for each aircraft type and govern safety of flight. In 
order to carry enough fuel to complete the rriission, and have the requirecTreserve 
fuel at landing, there may need to be a reduction in the number oT AXBT's or 
other sensors to be deployed in order to remain within NATOPS limitations. 

• Aircraft airspeed and altitude. These factors affect aircraft fuel consumption. 
Aircraft engines operate more efficiently at high altitude yet the ADDAS requires 
a lower altitude Tor effective operation. A compromise is made. Normally, 
transit to and return from the survey region will is accomplished at high altitude 
(> 10,000 ft MSL) and at an optimum speed ('~310 Kts indicated airspeed). 
Once on-station, the aircraft will be operated at loiter airspeed (ISO to 210 kts) 
and altitude will normallv be between oOO and 1000 ft AGL (above ground levelT 
This will cover as larje an area as possible, and provide for optirnum ADDAS 
performance. .MrcratT altitude affects the time it takes for the AXBT to impact 
the water after deployment from the aircraft and the line-of-sight radio 
transmission range. At 'higher altitudes, the transmission range is increased and 
interference from' other sensors on the same RF channel may affect the AXBT 
transmission. Airspeed and altitude must be considered in fuel planning i, weight 
and balance)^ for the mission. .Additionally, safety of flight requirements as 
specified bv NATOPS and squadron SOP (standard operating procedures) may 
require higher altitudes for certain evolutions such as pilot changes. 

• .AXBT fall duration through the water, for shallow ('300 m» and deep (~60 m) 
AXBT's. 200 and 500 seconds, respectively (Sessions, er ai.. 1976, Bane and 
Sessions, 1984). 

• Nurnber of AXBT channels. Only three channels, 12. 14 , and 16 (corresponding 
to 170.5. 172.0, and 173.5 Vihz. respectiveiv. ) are currently manutactured (Gent, 
1982). Two. or more AXBT's .of the same'channel must not oe deploved vithin 
the ;ine-of-sight RF transmission range. The conflicting signals can trigger the 
-VI K9 and result in the loss of a profile. It is basicaifv a iunction oGaircraft 
aititude and speed. 

• External radio interference on .AXBT carrier channels. Once the .MK9 is 
initialized and aw^aiting the AXBT signal, spurious radio transmissions on the 
receiving channel can pre-trigger the .VlK9 and the digitized profile may be lost. 



D. SIMILAR SYSTEMS IN USE TODAY 

In addition to the Scripps system, the Naval Ocean Research and Development 
Activity (NORDA) also has a similar system, the Airborne Data Acquisition and 
Processing System (ADAPS) (Miles, 1984). ADAPS is currently used by 
NAVOCEANO for work in acoustics. These systems were designed for specific 
missions and are multi-channel capable. The organizations have the advantage of a 
technical staff and maintenance support, and as a result these systems are 
one-of-a-kind in design and operation. The NPS ADDAS, in comparison, utilizes 
totally off-the-shelf hardware. It also has the added capability of being employed in 
shipboard XBT surveys as well with a modified software package and different 
digitizing circuit board for the MK9. However, the ADDAS is limited in that it is only 
single-channel capable at the present. With further advances in expendable sensors, 
off-the-shelf digitizing units, and computer systems capable of multi-tasking, systems 
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such as the ADDAS, ADAPS, and the Scripps system can become available to the 
Naval Oceanographer and research community to more fully describe the oceanic 
environment. 



16 



III. OCEANOGRAPHIC APPLICATION OF THE ADDAS IN THE 

OPTOMA PROGRAM 



Near-real-time synoptic mappings of the mesoscale eddy field are a key facet in 
the observational network of an ODPS. These mappings, using data acquired through 
the OPTOMA Program's AX3T survey flights, provide nowcasts of the mesoscaie 
variablity in the OPTOMA domains which are used for shipboard survey planning, 
dynamic modei initialization, model forecast verification, and assessment of the 
interannual variability in the domains. 

The NOCAL (Northern California) and CENC.AL (Central California) 
subdomains ire the major survey regions for the OPTOMA program (Figure 5.1). 
Eacn subdomain is an open ocean 'egion m Tie California Current System iCCS) 
about 150-km square with water depths of die order of -i km. The \OC.\L 
subdomain, 150-km downstream from the Mendocino Escarpment and about 150-km 
offshore, is the primaip' survey subdomain. Its bottom topography is characterized by 
a gently siooing seatloor and no seamounts, although seamounts, the Mendocino 
Escarpment, and the continental margin are nearoy. NOCAL is the simplest domain 
that can be found off California for a region this size (Mooers, 1986). AXBT survey 
flights generally expand the subdomains to about 260-km square in order to acquire 
sufficient data to initialize the boundary conditions for model runs. 

A. AXBT SURVEY FLIGHT PLANNING 

As with any oceanographic survey, planning is essential for the successful 
completion of the mission. Formal tasking for research flights is requested through the 
chain of command from the research sponsor. AXBT survey flights are then scheduled 
on a not-to-interfere basis and are contingent upon the operational tasking of the 
patrol wing. Once flight time has been authorized, coordination with the Patrol Wing 
Schedules Officer, preferably two to three months in advance of the date of the 
proposed survey flight, should begin. Flexibility is the key, and periodic contact should 
be maintained at least bi-monthly as changes in operational commitments (within both 
parties) may require the rescheduling of the flight. One month prior to the flight, a 
draft flight plan is prepared, and required mission personnel, supplies, equipment, and 
transportation are identified. Within two weeks of the flight, the Schedules Officer 
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Figure 3.1 OPTOMA Program survey region subdomains. 
Isobaths are in meters. 
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should know which squadron will be tasked to provide the services. Coordination is 
now primarily with the squadron. 

The proposed flight plan should be discussed with the squadron Operations 
Officer and Schedules Officer. Changes for safety-of-flight, number of AXBT's to be 
deployed, and the planned flight itinerary should be included in the revised flight plan. 
Additionally, squadron personnel should be requested for unloading and preparing the 
AXBT's for the flight and loading them onto the aircraft. 

B. AXBT SURVEY FLIGHT EXECUTION 

OPTOMA survey flights have generally deployed between 70 and 96 AXBT's per 
mission. Once on station in the survey subdomain, AXBT's are normally launched at 
15 nautical mile (NM) intervals, alternating shallow and deep AXBT's. The SSQ-36A 
AXBT is manufactured in two factor-'-set depths. 305 m (currently the Na\w standard 
AXBT) and 760 m (primarily used for research). The shallow AXBT's are digitized in 
real-time and stored by the ADDAS, while the deep AXBT's are recorded on tape for 
post-flight digitizing. A flight log is maintained with specific information about each 
drop: time, latitude, longitude, remarks. All shallow AXBT's are of one channel, while 
all the deep AXBT's are of a different channel. This limits the radio interference 
between .AXBT drops since only two channels. 14 and 16, can be used in the 
OPTOMA survey regions. (Channel 12 is not used as it is the California Forest 
Service emergency broadcast channel.) 

With the survey completed, the stored AXBT profiles are read and the desired 
output variables; e.g., sea surface temperature (SST), mixed layer depth (MLD), depth 
of the 8°C isotherm (DSC), sound speed at a chosen depth (CSZ), or dynamic height 
relative to a specific depth (DHZ), are selected. These variables are then entered into a 
statistical objective analysis routine based on an application by Carter and Robinson 
(1981). A discussion of the objective analysis technique as it is applied in the ADDAS 
software is presented in Chapter IV. This technique uses a two-dimensional correlation 
function based on the statistics of the data set. The five most correlated values within 
a 55 km radius of an interpolation point are then used to compute a value at the 
interpolation point. Contoured maps of the selected fields are then generated and 
printed on the plotter. These plots are available prior to landing so that an immediate 
debrief can be accomplished. Upon return to NFS, the recorded deep AXBT's are 
digitized and added to the data set. The objectively analyzed fields are then generated 
on the IBM 3033 mainframe using both the shallow and deep AXBT's. 
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The most recent AXBT survey flights were OPTOMA 18P consisting of two 
consecutive flights on 31 October and 2 November 1985, and OPTOMA20P on 16 
March 1986. Summaries of each mission and the data acquired are presented below. 

1. OPTOMA 18P 

Aircraft and flight crews were provided by Patrol Squadron Forty-Six, 
assigned to Commander, Patrol Wing Ten, for the 31 October flight and Patrol 
Squadron Ninety-One, assigned to Commander, Reserve Patrol Wing Pacific Fleet, for 
the 2 November fight. Both are based at NAS Moffett Field. The flight tracks were 
executed as planned (Figure 3.2). A total of ISO AXBT's were deployed on the two 
flights. Spacing was at 15 NM inrervals, alternating deep (Ch. 16) and shallow (Ch. 
14) AXBT's. On-station ground speed averaged 200 kts for each flight and the altitude 
varied between 500 and 800 ft AGL, depending on the low-level visibility. The logged 
positions were obtained from the aircraft inertial navigation system with an accuracy of 
: NM km). 

u. Mission Critique 

Two different models of the P-3 aircraft were utilized on the fights: a P-3C 
on 31 October, and a P-3B on 2 November. The following differences are noted 
between the two models for ADDAS installation: 

• The 28VDC power connections in the P-3C were made to terminal board 426 
(TB426) located lust aft of sensor station 2 (SS2) bulkhead. The terminals used 
were CIO ( + 28vDC power SS2) and Cll (ground). In. the P-3B, power 
connections were made to TB431 located above SS2. Terminal A24. provided 
+ 28VDC power while the ground was a stud located on the far right side of the 
terminal board. The sonobuoy pin connections were the same for both models of 
airciaft. (Please see Colton and Mooers, 1985, for a complete description of 
ADDAS installation.) 

• The P-3C model does not have a LORAN system incorporated into its avionics 
package. Thus, no LORAN antenna was available to be connected to the 
ADDAS LORAN unit, so all AXBT positions were obtained from the inertial 
navigation system and checked against tne O.MEGA radio navigation unit. 

• The physical internal layout of the P-3C is different from the P-3B. .There is 
much less room for the ADDAS and operator in the P-3C at sensor station 2 due 
to the bulkhead separating stations. The ADDAS operator had to work literally 
"shoulder to shoulder" with the sensor station 1 operator. 

The 31 October flight with Patrol Squadron Forty-Six went very well. 
Since the ADDAS LORAN was not used, the positions of each AXBT drop were 
taken from the inertial navigation and entered manually into the AXBT program. 

The 2 November flight proceeded very smoothly due to the experience 
gained by Patrol Squadron Ninety-One from previous OPTOMA flights. ADDAS 
installation was completed in a record time of forty minutes. AXBT positions were 
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Figure 3.2 Flight Plan, OPTOMA ISP. 
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obtained as on 31 October due to problems encountered on previous flights with RF 
interference from aircraft radios affecting the LORAN signal. The only notable 
difficulty on this flight was that a low-level (500 ft) stratus deck required an on station 
altitude of 800 ft AGL for flight safety. This increased the transmission range of the 
AXBTs. On the first southbound leg, a tailwind increased aircraft groundspeed to 235 
kts and resulted in pre-triggering of the MR9 due to overlapping AXBT signals on the 
same channel. On subsequent downwind legs, airspeed was adjusted to maintain 200 
kts groundspeed. Table 1 lists a summary' of AXBT expenditures for the two flights. 



TABLE 1 

AXBT EXPENDITURE SUM.MARY OPTOMA18P 



Flight 


.Area 


Cliannei 


Depioj'ed 


Failed ^ 


Recovered 


31 Oct 


CENCAL 


16 




6 


39 






5 /I 


Sd 


1 


40 


2 Nov 


NOCAL 


16 


-i6 


6 


40 






i. -t 


-5 


4 


41 


Two -day 


Total: 






lo 


160 



Notes; 

' Due to no radio signal transmission and late light-otf of the probe. 



b. AXBT Survey Results 

The uncorrected positions of the shallow AXBT's, from the 31 October and 
the 2 November flights, and the analyses of D8C, SST, and CS150 that were generated 
onboard the aircraft on the return leg from the 2 November flight were examined 
(Figures 3.3a, 3.3b, 3.3c, 3.3d, respectively). A compiled version of the objective 
analysis program, using a floating-point accelerator was applied. Only the shallow 
AXBT data (30 NM spacing) were combined to produce the analyses. The data were 
used as received and digitized by the ADDAS without editing or filtering performed on 
profiles in order to provide a real-time analysis on the return flight to NAS Moffett 
Field. Some quality checks are made by the ADDAS software. A complete summary 
of the acquired AXBT data is presented by Wittmann et al. (1986). 
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Figure 3.3a Shallow AXBT Drop Locations, OPTOMA18P. 
Origin 38"N, 126“W. Tik marks = 25 km. 
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Figure 3.3b Depth of the 8“C Isotherm (m), as in Figure 3.3a, 
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Figure 3.3c Sea Surface Temperature ("C), as in Figure 3.3a. 
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Figure 3.3d Sound Speed at 150 m (m/s), as in Figure 3.3a. 
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The inflight digitized channel 14 AXBT data were later combined with 
post-flight digitized channel 16 AXBT at NFS. These data were stored and processed 
on the IBM 3033 mainframe computer. Profiles were edited to remove spikes, and bad 
profiles were eliminated. Post-flight, the recovered AXBT locations were examined 
(Figure 3.4a). 

• DSC, meters (Figure 3.4b): This analysis presents a "topographic" map of the 
depth of the S^C surface, which is generallv accepted as a aooa indication of the 
location of the permanent thermocfine ancf pvcnocline for this region. From the 
DSC analyses, the streamfunction for the surface geostrophic flow' can be inferred 
with good accuracy. The major features were two cyclones and two anticyclones 
(indicated bv "C's"' and "A's", respectively), alternating along the coast. The 
cyclone at the far nortn and the anticyclone at the far south were fairly weax 
compared to the cyclone close to the coast near Point Reyes (diameter about 100 
km, depth variatio'n of over 60 m) and the stronger anticyclone pfl' Point Arena^ 
(diameter about 150 km, depth variation of over 60 m). The juxtaposition of 
these eddies contributed to alternating onshore and offshore flow. 

• SST, °C (Figure 3.4c): The main features are two cold core centers ( indicated by 
"K's”) near tire coast and a zeneraily weak positive orfshore temperature aradienc. 
The cold core center in the s'outh of the sampled reaion was elohaated, parallel to 
the coast, whereas 'hat m hie north of the 'sampled reafon (near Cape 
Mendocino) was more circular, slightly elongated oflfshofe. 

• CS150. m/s (Fiaure 3.4d): Sound speed was computed using a constant .salinitv 
of 33.5 ppt. The sound speed is fairiv constant with only a 6 m/s variaouity, 
although there is an obvious, correlation between sound speed patterns and the 
eddy patterns of the depth of tiie 8 dearee isotherm, ihe strong reature evident 
in FA'ure 3.3d mst to the south of the'northernmost cold center was a prooiern 
with fne data which was eliminated in the post-flight data processing. 

2. CPTOM.420P 

A combined flight crew from the Naval Air Reserve Center and Patrol 
Squadron Ninety-One at NAS Moffett Field performed the mission. The flight plan 
was executed as shown in Figure 3.5. The ten flight tracks were spaced 15 NM apart. 
Ten AXBT's were deployed along each track at a 15 NM spacing, with the exception 
of only eight AXBT's deployed on tracks 1 and 10 (no AXBT's were deployed at the 
comers of the survey region), for a total of 96 AXBT's. Groundspeed averaged 220 kts 
and on-station altitude varied between 500 and 2000 ft AGL for safety of flight due to 
turbulence and scattered convective activity and thunderstorms in the area. 

AXBT's were to be deployed alternating channel 16 (deep, 760 m probe) and 
channel 14 (shallow, 305 m probe). Due to inventory error, all AXBT's for this flight, 
both channel 14 and 16, were shallow (305 m). This was not realized until the first few 
channel 16 AXBT's were deployed and their time of fall observed. AXBT's were then 
launched at approximately four minute intervals, alternating between channel 16 and 
14 AXBT's. AXBT positions were obtained from the aircraft inertial navigation system 
(2 NM accuracy) and entered in the flight data log. The positions were not corrected 
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Figure 3.4a Recovered AXBT Drop Locations, OPTOMA ISP. 
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Figure 3.4b Depth of the 8° C Isotherm (m), OPTOMA18P. 
Contour Interval 10 m. 
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Figure 3.4c Sea Surface Temperature TC), 0PT0.MA18P. 
Contour Interval, 0.25 C. 
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Figure 3.4d Sound Speed at 150 m (m/s), OPTOMA 18P. 
Contour Interval, 0.5 m/s. 
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Figure 3.5 Flight Plan, OPTOMA20P. 
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for navigation error. The channel 14 positions were manually entered in the AXBT 
acquisition program and the signal was digitized in real time by the ADDAS and 
stored on diskette. Both channel 14 and 16 AXBT analog audio signals were recorded 
on magnetic tape using a stereo hi-fi VCR. Channel 16 AXBT's were digitized on 
return to NFS. 

c Mission Critique 

The flight proceeded very well even though the weather was less than 
optimum. The survey region was dominated by low stratus and cumulus clouds as well 
as a few towering cumulus. Moderate turbulence was experienced throughout the 
duration of the flight. Amditionaliy, scattered rain squalls required deviation from the 
planned flight tracks for safety of flight. Winds at flight level were southeastward at 30 
knots, and the sea state was 3. The navigation team did an excellent job and the 
resulting .AXBT droos ended up ver>' nicely spaced and close to planned drop positions. 

.A malfunctioning channel 16 sonobuoy receiver on the aircraft distorted 
the .AXBT audio earner signal The ADDAS oscilloscope provided a visual reference 
for the Inilight Technician (IFT) to troubleshoot and repair the receiver. Fortunately, 
only 8 channel 16 .AXBT audio signals were lost and could not be digitized by the 
Sippican Mark 9 unit upon return to NFS. Due to a minor problem witii :lie 
lofargram recorder, only 6 of the 8 AXBT traces were obtained and these profiles were 
hand digitized. 

The use of a stereo hi-fi video cassette recorder (VCR) in the audio mode 
proved to be a worthwhile addition to the ADDAS system. With its capability of 
recording up to six hours of audio on a T-120 VHS cassette, the operators only had to 
change the tape once as compared to previous flights where the AXBT audio signal 
was stored on standard audio cassettes with only 45 minutes of recording on each side 
of the tape. No profiles were lost during tape changes or because a tape ran out. The 
VCR was mounted in the ADDAS cabinet on a layer of 3/4-inch foam rubber, which 
helped insulate it from vibration caused by turbulence encountered during the flight. 
No degradation in the recorded signal was noted. 

The internal ADDAS inverter is rated at 5.6 amps. With the inclusion of a 
VCR in the system, the inverter was running at 6.3 amps continuous. This did not 
appear to be a problem and the inverter did not overheat or malfunction. To be on 
the safe side though, once data acquisition was completed, the VCR was disconnected 
from the circuit when the plotter was to be used. 
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The Sippican MK9 digitizes the AXBT audio signal ten times each second. 
MK9 output was sampled for three-tenths of a second, producing three data points. 
The median value of the three digitized points was stored. This step was included in an 
attempt to eliminate non-representative spikes in the temperature trace caused by RF 
or other interference in the signal. This feature appeared to work well, with the 
exception of long duration breaks in the signal, possibly caused by sea wash-over of 
the AXBT transmitter. A minimum of post-flight editing was required on the profiles. 
(A compilation of the AXBT data is presented by Ciandro et al., 1986.) Table _ lists a 
summary' of AXBT expenditures for this mis.sion. 



TABLE 2 

AXBT EXPENDITURE SUMMARY OPTOMA29P 



Channel Depioved F aiied^ 


Lost" 


Deleted'* 


Recovered 


M X8 1 


0 


0 


47 


16 48 6 


- 


1 


39 


Total: 96 


•> 


1 




Notes: 








^ Due to no radio signal transmission 


and late light-off of 


the probe. 



^ Profiles lost due to absence of lofargram trace from which to hand digitize. 
^ Bad profiles deleted during post-flight editing. 



b. AXBT Survey Results 

A representative sample of AXBT traces show the variability of the thermal 
structure in the survey region (Figures 3.6a and 3.6b). The AXBT drop positions, and 
the analyzed fields of DSC, SST, and MLD were produced by the ADDAS on the return 
leg of the flight (Figures 3.7a, 3.7b, 3.7c, 3.7d). Only channel 14 AXBT data (30 NM 
spacing), as received and digitized by the ADDAS, were used. The fields were 
generated by a compiled version of the objective analysis program using a 
floating-point accelerator. 

Post-flight digitized channel 16 AXBT data were combined with the inflight 
digitized channel 14 AXBT data (Figure 3.8a). These data were stored and processed 
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Figure 3.6a Representative Temperature Profiles, OPTOMA20P. 
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Figure 3.6b 



As in Figure 3.6a. 
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Figure 3.7a Shallow AXBT Drop Locations, OPTOMA20P. 
Origin 38°N, 125"W. Tick marks = 25 km. 
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Figure 3.7b Depth of the 8®C Isotherm (m), as in Figure 3.7a. 
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Figure 3.7c Sea Surface Temperature (°C), as in Figure 3.7a. 
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Figure 3.7d Mixed Layer Depth (m), as in Figure 3.7a. 
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Figure 3.8a Recovered AXBT Drop Locations, OPTOMA20P. 
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on the IBM 3033 mainframe computer at NFS. Remaining spikes were removed and 
bad profiles were deleted from the data base. 

• DSC. meters (Figure 3.8b): The most notable features are the major cvclonic 
(counterclockwise, indicated by a "C" ) center in the northwest with its associated 
flow zone extending to the south-southeast, as indicated bv the shoaler depths, 
and the very' large anticyclonic (clockwise, indicated "by an "A") center 
characterized 'by the depre'ssion in the southeast. Additionallv. half of two 
smaller cyclonic flows are observed in the north as well as a small cyclone near 
the center of the region. These features indicate an onshore jet off Ft. Arena. 

• SST, ”0 (Fiaure 3.8c): Major features of interest are the cool zone in Che 

northwest and the warm center (indicated by a ’'W"). correspondina to che larae 
cvcione in- che northwest and the large anticyclone to the southeast) respectiveiv.. 
There is also evidence of a cool center ("K")m the north corresponding to one of 
the two smaller cyclonic flows. 

• MLD. meters (Fiaure 3.Sd): Of interest in this analysis is che wide range of 
values for MLD. ifom less than 16 m to more than 68 'm. The mean .VILDTs iO 
m. This variabilitv could be che result of the larae atmosohenc occluded 
low-pressure svsterh chat dominated the reaion on the previous dav. As the 
system passed.' tne winds shifted, from generally northeastward to southeastward 
a't 30 to -0 knots and the surface layer was apparently in a state of rapid change 
durina the survev period. .All AXB i profiles 'showed a ven.' isothermal surface 
la 3 /er.' .here was some association of the patterns in MLD' with chose oi' D8C; 
fo'r example, me MLD had prominent minima associated with both the cyclone 
'n the northeast and the anticyclone in the southeast. 
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Figure 3.8b Depth of the 8°C Isotherm (m), OPTOMA20P. 
Contour interval, 5 m. 
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Figure 3. Sc Sea Surface Temperature (!C), OPTOMA20P. 
Contour interval, 0.1®C. 
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Figure 3.8d Mixed Layer Depth (m), OPTOMA20P. 
Contour interval, 2 m. 
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IV. DEVELOPMENT OF AN AIRBORNE OCEAN SAMPLING 

STRATEGY 



There are two distinct considerations that are related to the successful areal 
mapping of oceanic variables; 1) The sample density in the surv'ey region; and 2) The 
application of the objective analysis technique to map these sam^ples. 

A. SAMPLING CONSIDERATIONS 

Given an ocean region to be synoptically mapped in some variable, and the 
aircraft asset(s) and air-depioyabie sensors to measure this variable, what airborne 
sampling strategy should be employed to acquire a representative sample of the field 
that will optimize the information returned with the minimum amount of time and 
euort expended to accurately map this variable/ 

The ocean is a vast, temporally and spatially varvang mhomogeneous domain. 
Following a simiiar development by Kelley (1976). tliis domain can be described by 
basic probability theoip/ as a sample space. S ISj, Sn,...5^), of n observations where 
-S; = (x:, vj, t;}, i= l....n, is an dement of the sample space deiined by latitude, 
longitude, and time. Each element can be further described by an array, V; : 
(Kj j(z), K- 2(z),...Pj j(z)), j=l,...,m, of specific oceanic variables (biological, chemical, 
dynamic, or physical) varying with depth (z). Each F- ;(z) can be defined as a 
continuous random variable which can have more than one potential value as the 
chance result of the sampling and that it is assumed to be representative of the oceanic 
variable being sampled. Due to the inherent capabilities of sampling instruments and 
systems such as the ADDAS, these oceanic variables are sampled at a specific depth 
increment and can be represented as discrete random variables. This provides a 
framework for the development of an "optimum" sampling strategy that can be applied 
in survey operations. 

Deciding which variable to sample is a rather simple task due to the availability 
of air-deployable sensors. The measurement of ocean temperature is most practical 
through the use of the SSQ-36 AXBT as it is a proven sensor in the Navy inventory. 
Synoptic mappings of the temperature field at selected depths can be converted to a 
sound speed field, using a climatological or representative salinity profile, for use in 
acoustic forecasting for ASW operations. Soon-to-be-operational AXCTD's and 
currently available AXSV's could also be used to map the sound speed field . 
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Having chosen the variable to be sampled, the "external" considerations of the 
spatial and temporal scales of the variable, and the "internal" considerations relating to 
aircraft, aircrew, and sensor limitations will be addressed. 

1. External Considerations 

a. Spatial Scales 

The ocean contains a multitude of spatial scales ranging from centimeters 
to thousands of kilometers. In sampling any oceanic variable, prior knowledge of the 
scales that are "significant" to a specific mission are required. Previous operations in 
the survey region, satellite imager\% and/or climatology may be available. Additionally, 
a knowledge of the dynamical properties of the field must also be considered. 

The Rossby radius of deformation, R, is a length scale of fundamental 
importance in both atmospheric and oceanic dynamics for the behavior of rotating 
density-strathied fluids subject to gravitational restoring forces. For a homogeneous, 
shallow layer of iluid, it is the distance over which the gravitational tendency to render 
the free surface Hat is balanced by the tendency of the Coriolis acceleration to deform 
the surface (Pedlosky, 1979; p, ~8), It is the pertinent length scale for barotropic 
phenomena (Emery, ei al.. 1984). Since the ocean is a continuously stratified Iluid. 
normal mode theory can be applied to a stratified layer and each mode can be 
represented by a Rossby radius. The barotropic or external Rossby radius (R) was 
defined above and is a function of only the water depth and latitude. The internal or 
rtth mode baroclinic Rossby radius (Ri) depends on the vertical stratification and is a 
natural scale in the ocean often associated with fronts and eddies (Gill, 1982; p. 207). 
For mid-latitude oceans it is the order of 10 to 30 km. The horizontal dimensions of 
such quasi-geostrophic features should scale as the internal Rossby radius, where the 
scale is defined as the wavelength over 2n (Emery, et ai, 1984). 

b. Temporal Scales 

Based on the spatial scales of the quasi-geostrophic features, the temporal 
scales of variability are not a factor in "synoptic" mappings. With phase speeds the 
order of 2 to 5 km/day, the distance that the features will have propagated over the 
duration of the survey flight is small compared to the sampling grid cell length. 
However, temporal scales may have to be considered if the survey domain is too large 
to be covered in a single flight. 
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c. Aliasing 

The well known problem of aliasing occurring in time series analysis also 
occurs in spatial sampling. For a one-dimensional' sampling, the Nyquist theorem 
states that some function can be reconstructed from samples taken not more than 
one-half wavelength apart. This can be applied to a two-dimensional field as well. 
Identify the highest frequency (shortest wavelength) of significant variability and 
sample at twice that frequency (one-half of that wavelength). If this is not feasible and 
a larger sampling distance is chosen, aliasing may be a problem from scales that are 
smaller than the chosen sampling wavelength and may account for a percentage of the 
noise in a sample value. 

2. Internal Considerations 

With the spatial scale of the features determined, and thus the maximum 
acceptable distance between observations to adequately represent the significant scales, 
the capabilities of the sampling platform and its sensors must be taken into account. 
These were presented as extrinsic limitations to the ADDAS in Chapter 2, and are 
applicable here. They are summarized: 

a. Aircraft Limitations 

» Weisht and balance. Aircraft are constrained bv a maximum gross weight for 
take'otf as prescribed bv the respective NA i OPS '^lanuai. There may have to be 
a trade-off in the number of sensors to be deployed to remain within NATOPS 
limitations. 

• Airspeed and altitude. These directly affect fuel management. Airspeed and 
altitude will be governed by safety-of-flight and ADDAS/sensor limitations. 

Additionally, physiological and safety factors must be considered. All Navy flight 

operations are governed by the OPNAVINST 3710.7 series (NATOPS General Flight 

and Operating Instructions). This specifies limits on mission endurance and crew rest 

policies. It may be further modified by individual squadron standard operating 

procedure (SOP). These will be the limiting factors in mission duration. 

b. Sensor Limitations 

• AXBT deployment time from the aircraft, a function of altitude. 

• AXBT fall duration through the water, 200 and 500 seconds for deep and 
shallow AXBT's, respectively. 

• Three AXBT channels available. Two AXBT's must not be deployed within 
line-of-sight transmission range of each other, a function of altitude. 
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B. GENERAL APPLICATION 

The thought process followed in developing an "optimum" sampling strategy is 
presented in flow diagram form (Figure 4.1). Each step is contingent on the previous 
step. Maximum observation spacing is generally given as a requirement to gain an 
adequate representation of the oceanic process. Aircraft maximum gross weight is also 
a limiting factor. Trade-offs can easily be made between the size of the survey region 
and the maximum number of AXBT's that can be carried. If the size of the survey 
region is then too large, two or more consecutive survey flights may be completed as 
long as the temporal scales of the region are considered. One option, is to sample the 
larger do.nain twice. This can be accompiisiied by sampling in a "cneckerboard" 
pattern, with sample spacing no more than twice the maximum observation spacing, 
with the first flight deploying AXBT's in the "red" squares, and the second flight 
deploying in the 'blacK." squares. Another option is to complete two or mere surveys 
at a smailer-than-maximum spacing in adjacent subregions of the sample domain. One 
drawback to sampling with two or more tliguts is that the real-time aspect of the 
synoptic mapping will be lost. This may be acceptable, however, depenumg on the 
requirements directing the survev. 

C. TiiE OBJECTIVE ANALYSIS TECHNIQUE 

The second area of consideration in developing an optimum sampling strategy is 
the proper application of the objective analysis technique. This will result in fields that 
are representative of the data acquired, and available in near-real-time with systems 
such as the ADDAS. 

Based on the Gauss-Markov theorem from statistical estimation and applied to 
atmospheric data by Gandin (1965), objective analysis was developed by Bretherton, et 
al., (1976), with application to oceanic data. A more generalized implementation, as 
adduced by Carter and Robinson, (1981), is applied in the OPTOMA Program. 

1. Theory 

The basic assumption in objective analysis is that the field from which the 
observations were acquired has zero mean and a known covariance function. 
Following Bretherton, et al, (1976), and incorporating both a space and time 
variability, 6p is the value of a scalar variable 0(;r,y', r), at a general position p = 
(x,y,t), from measurements (pj at N observations pj ( i = 1, . . ., iV). Each 
measurement, (p-, is presumed to consist of some true value, 0j, and some random 
noise, Cj: 
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Figure 4.1 Optimum Sampling Strategy Flow Diagram. 
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(pi = 0(Pi) + Ci, (i = 

where the errors, £•, are uncorrelated with each other and with the true value, 9-, but 
have a known variance, E: 



Ci6i - 0, 






- -ESij. 



for ( i, j = 



and the overbar indicates an expected value (or average over a large number of 
realizations). The errors may be sampling error or instrumental error. Systematic or 
calibration errors are not permitted. 

The least squares optimum linear estimator of 9j is 



A 

8. 



N 

y c. 



iV 



'PI ^ - ^ii 
1=1 i=i 



«ij 



wher: 



Ajj = ^ p[- ?j ) + 

is the covariance matrix between all pairs of the sampled data, and A'* is 'he inverse of 
A. and 



Si " Vi " 

is the covariance between the estimated value 9p and the \th observation. The variance 
of the error in 9p is 

A ^ 

( Sp - 0p )' = Cpp ■ S Cpi Cpj Aij>. 
i,j= 1 

which requires only the location of the data points, a noise level E, and some 
correlation function F{R). 

Carter and Robinson, (1981), present the idea of applying the objective 
analysis equations at each interpolation point. Four specific parameters are identified 
that determine the selection of observations that will affect an interpolation. These are; 

• Maximum time increment between observations, 

• Radius of influential points, R = [ ( Ax - c^^Ar + A>’^ ] , 

• A dominant phase speed, c^^, 

• The maximum number of influential points. 
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With application to synoptic mapping, only the radius of influence and number of 
influential points are used. It is best to use observation points near the interpolation 
point.- By specifying a maximum range at which observations can effect an 
interpolation, this will eliminate the influence from distant points. Additionally, the 
specified maximum number of points limits only those observations within the radius 
that are the most highly correlated (closest to) an the interpolation point. The time 
increment and phase speed are disregarded because the propagation distance of the 
quasi-geostrophic features (phase speeds of the order 2 to 5 km/day), is negligible over 
the duration of a survey flight (~8 hrs). The radius of influence is now only the 

O 1/ 

distance between the interpolation point and the observation /? = ( A.r" + AjT- 
2. ADDAS Software Application 
a. Input Parameters 

Having acquired the observations, the following parameters must be 
specified before performing the objective analysis: 

• Some form of the correlation function. F( /? ), 

• The percentage of the observed variance attributed to noise, 

• The maximum radius of influence. 

» The limiting numoer of iniluentiai observations. 

• The objective analysis interpolation points, = { ^n’-^n ~ ). 

• A maximum autocorrelation value, 

These parameters may be incorporated into the program as constants based on 
previous experience or determined for each observation field. A flow diagram (Figure 
4.2) summarizes the decision process. These parameters are presented below. 

(1) Correlation Function and Noise Variance Level. 

There are two distinct uses for the correlation function in the objective 
analysis technique: 1) to determine the correlation between an interpolation point and 
an observation, and 2) to compute the elements of the autocorrelation matrix between 
all pairs of observations selected to influence an interpolated value. 

The correlation function can be defined in any number of ways 
depending on the user's requirements. For a space-time correlation of asynoptic data, 
as is used in the statistical model, the form of the correlation function is 

F(i?2) = Y( 1 -a/?2)exp(-pf?2), (4.1) 

where 

R2 = ( Nx - c^At )^ H- ( Ay - CyAt )^, the influential radius, squared. 



52 




Figure 4.2 Flow Diagram of Objective Analysis Inputs. 
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c^, Cy: X and y phase velocity components, respectively, 

Ax, Ay: spatial difference between interpolation and observation points. 

At: time difference between observation and analysis tim^es, 

a, P, j: parameters computed bv an iterative nonlinear least squares fit to 
the function. ' ^ 

A simple, yet effective form is employed in the ADDAS objective 
analysis software for synoptic mappings. 

F( R2 = 1 1 - R.2 ^:qxv( -R2!2 ), (4.2) 

where 

R2 = Ax — ^ Ay" ) (Rl^.^y.p', the influential radius, so uared. 

•\n advantage in using equation ^.1 in synoptic mapping 'Jot and c„ 
both equal zero) is that an estimate of the noise variance in the observations can be 
made if not already known. First, a bilinear trend is removed (if desiredl from each 
observation. It is of the :brm <0: - ,\x; - 3v’; -C. vvmere A. 3. and C are either 
specified or computed by a least squares method. The mean, .ii, and the variance, (?", 
of the observations are computed next, and the mean is removed from each 
observation to satisfy the basic assumptions. 

The distances between each observation and every other observation, 

Dij = + (ypyj)2]'/^(i,j = 

are computed and binned as if isotropic. Bin size should be greater than the minimum 
inter-observation distance ( D-j, i j ) in order to accumulate a statistically significant 
number of values in each bin. The first bin should be from 0 to 1 km, to contain those 
D;; when i = i, and other closely-spaced observations ( D;; ^ 1 km ). The number of 
bins, nbin, should cover all the inter-observation distance scales. As each Djj is located 
in its respective bin, b, the total number of elements in that bin, NB^, is incremented 
by 1. The autocovariance is computed and summed for that bin ( ^ <pj<pj )• 

Once all the pairs of observations have been binned, the final autocovariance for each 
bin, is computed by normalizing each by the quantity 

(C|, = C0|,/(WS|,ff-)). 
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The midpoints of the bins, BM^, for b = 2,...,nbin, are computed. 
These values, along with the total number of observations ( NB^ ), the normalized 
correlations ( ), the variance (<y^), and the limiting number of observations { N/2), 

are applied in an iterative nonlinear least squares fit to the correlation function 
(Equation 4.1). The values of a, p, y, and the zero-crossing distance are computed. 
Some quality checks are made on the data. If the total number of values in either bin 
2, 3, or 4 is less than N/2, there are not a statistically significant number of 
observations used in computing the correlation, and the least squares fit is not 
attempted. Additionally, the application of the least squares routine .Hor this .Horm of 
the correlation function requires at least three points to fit a curve, and these points 
must also occur before the second zero-crossing of the correlation values. The 95% 
confidence limits, based on a test for white noise and assuming a normal distribution, 
can be computed for each bin using the formula Ci^ ± 1.96/( NBi^^ '^'y (Jenkins and 
Watts, 1968: pp. 187-1881. 

If a successful fit is achieved, the signal and noise parts of the variance 
can be estimated by extrapolating the fitted correlation function to zero lag. Assuming 
perfect correlation at zero lag (C= 1.0) and that the variance contains some 
percentage of noise, if the extrapolated value. is less than l.O. the signal variance 

can be computed, 

and the noise variance is 

The noise variance must be specified as a fraction of the total variance for input into 
the objective analysis routine. 

(2) Maximum Radius of Influence. 

i?/max selected to be less than or equal to the zero-crossing distance 
of the correlation function. This ensures that only observations that are positively 
correlated with the interpolation point are used. 

(3) Maximum Number of Influential Observations. 

In surveys with a large number of observations, Ny^^^ limits the 
number of observations within the radius that will be used in an interpolation. These 
observations will be those that are the most highly correlated with (i.e. closest to) the 
interpolation point. 
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(4) Grid Selection. 

For synoptic mappings in the OPTOMA Program, the survey domains 
are rectangular in shape, with observations distributed regularly over the domains. The 
objective analysis interpolation points, V}/j^ — { }. ~ are organized in a 

rectangular pattern encompassing the domain, with equal spacing in both the local x 
and y directions. The grid spacing should be less than one-half of the zero-crossing 
distance of the correlation function to adequately resolve the scales of the sampled 
field. The interpolation points should lie within the dimensions of the domain in order 
to provide enough observations for interpolation at the boundaries. 

(5) Maximum Autocorrelation Value. 

/I max computing the autocorrelation matrix. A, at each 

interpolation point. Once the observations are selected that will influence an 
interpolation, the autocorrelation of each point with every other point, .A;j, is 
computed using the autocorrelation function F{ R ) previously determined. If for a 
given i, the absolute value of A;j, is greater than 

observation will be discarded. This removes observations that are too closely spaced to 
one another and helps maintain a diagonally-dominant matrix. 
b. The Objective Analysis 

Once the input parameters are specified, the following steps are performed 
for each interpolation point, 

• Selection of the influential observations, 

• Generation of the autocorrelation matrix, 

• Computation of the interpolated value and error limit. 

These are discussed below. 

(1) Selection of Influential Observations. 

The influential radius is computed between the interpolation point and 
each observation. If the radius is less than or equal to correlation 

between the observation and the interpolation point = F{ R), is computed. Once 
the observations have been selected, they are sorted in order of decreasing correlation 
(increasing distance from an interpolation point). Only the first of the sorted 

observations and their respective correlations will be used. Next, the autocorrelations, 
Ajj, between each of the selected observations are computed. Only the best correlated 
and distributed observations with the interpolation point remain. 
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(2) The Autocorrelation Matrix. 

The autocorrelations, A-, between the remaining selected observations 
are calculated using F{ R ), and become the elements of the autocorrelation matrix. 
The percentage of the variance due to noise is then added to each element of the main 
diagonal of A. This biases the the matrix towards diagonal dominance. The matrix is 
then inverted. This provides the weights to compute the interpolated value. The 
observation correlation vector between each selected observation and the 

interpolation point is computed using F{ R ). 

( 3 ) The interpolated Value and Interpolation Error. 

The vaiues of the inverted autocorrelation matrix. A'^ and the 
correlation vector, C^; are used to compute the interpolated value, B^, 



V .V 

8a- I <^ai I ^ 

i -1 1=1 ■ ■ 



and the RMS interpolation error is 
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The interpolated value must now be be brought back to its original 
scale. If a bi-linear trend was removed, it must be added back, as well as the mean. 
The final result, R^, is 

^n ^ ®n "^^n ®-^n 



which is stored, along with E for later contour plotting. The process is again 
executed for the remaining interpolation points. 

3. Possible Problem Areas 

From the numerous runs of the objective analysis that have been made in 
gaining an understanding of the technique and processing data, two interrelated factors 
have been found that determine the caliber of the objectively analyzed fields: data 
density and quality, and the noise variance estimate. 

Generally, the more robust the sample field, the better the quality of the 
resulting objective analysis. However, this may not always be the case. Too many 
selected observations too close to each other (high autocorrelation values) may result 
in an autocorrelation matrix that is not diagonally dominant, and an ill-conditioned 
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(singular or near- singular) matrix may result. If this happens to be the case, the 
interpolated value may be grossly non-representative of the data it was computed from, 
and the RMS interpolation error will be high. By specifying ^ 

certain value (limiting the minimum distance between observations), this will limit the 
magnitude of the off-diagonal elements and result in a properly conditioned matrix, 
which should produce acceptable results. 

The estimated noise variance is another factor that can affect an interpolation, 
although the erfect will not be as great as improperly specifying the maximum 
autocorrelation value. If the noise variance is specified to be too small, especially in 
sparse data fields vvith a high noise vanance, this may also cause the autocorrelation 
matrix to become ill-conditioned. As the noise is added to the main diagonal of the 
matrix, it also biases it toward diagonal dominance so a successful inversion can be 
accomplished. With a higher density of observations, smaller bin sizes can be used to 
compute the autocorrelation values and more points will be available for the least 
squares fit. This could result in a oetter estimate of that part of the variance 
attributabie to noise. 

.Another name for the objective analysis techniaue is optimal interpolation. 
With a knowledge of the statistics of the data and a correlation function, the fraction 
of the variance due to noise can be calculated. These factors, along with a specified 
maximum autocorrelation value and the observations will result in an optimal 
interpolation of the observations to each interpolation point, and an acceptable 
representation of the sampled field. Failure to properly specify these parameters may 
result in erroneous mappings due to the problems of non-singular matrices and 
observations that are too closely spaced. The effectiveness of a system such as the 
ADDAS for synoptic mappings is directly related to quality of the acquired data and 
the proper processing of this data. 

D. THE OPTOMA PROGRAM SAMPLING STRATEGY 

Applying the considerations discussed in the previous sections, and with the 
knowledge gained from earlier OPTOMA Program AXBT survey flights, the following 
key factors were considered in the development of the airborne sampling strategy in the 
NOCAL and CENCAL domains of the OPTOMA Program: 

• Size of region: The 250-km square covers the maximum area that can be 
surveyed by a Navy P^3 .aircraft based at Moffett Field Naval Air Station, 
California, on a single mission. 
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• Number of AXBT's: 96 AXBT's is the maximum number that can be carried by 
the P-3 to remain within specified gross weight limitations and still have adequate 
fuel, including reserve, to complete a mission in this size region and at this 
distance from base, and to provide an adequate resolution of the inesoscale 
processes. 

• AXBT spacing; This is a function of many variables. AXBT's are deployed at 
five minute time intervals ( ~ 28 km spacing). . This increment has been 
determined from previous flights to provide an optimal coverage of the region, 
and it takes into account the spatial scale of oceanic features to be observed, 
ADDAS response, AXBT fall time from the aircraft, AXBT transmission time, 
and AXBT-aircraft line-of-sight transmission range, for aircraft on-station 
groundspeed of 180 to 200 kts and altitude of 500 to TOOO ft AGL. 
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V. SIMULATED APPLICATION OF THE SAMPLING STRATEGY 



The sampling strategy was applied to "simulated" AXBT survey flights in 
"synthetic" oceanic domains. The oceanic variability in these domains was synthesized 
from the superposition of linear planetary Rossby waves (LPRW's). The simulated 
survey domains were 256-k.m by 256-km square, representing a nominal areal survey 
domain as in the OPTOMA Program. Synoptic mappings were generated from the 
simulated AXBT observations using the objective analysis technique. Some 
assumptions were made in applying the sampling strategy; these were: 

* No limit to aircraft endurance, airspeed, or payload capability, 

* Unlimited numoer of AXBT's, 

* No .AXBT failures. 

* No limitations on the ADDAS. 

* Temporal scales of the iieids were not considered. 

.Although not fully realistic, these assumptions provide a basis to exercise the objective 
analysis tecnnique on controlled data. 



A. DEVELOPING THE SYNTHETIC FIELDS 

A field of eddy-like features can be generated from the superposition of LPRW's. 
Starting with the non-linear, hydrostatic, Boussinesq equations, and using a 
non-dimensionalization and perturbation expansion in terms of the Rossby number, a 
quasi-geostrophic potential vorticity (QGPV) equation can be derived: 



where 
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(The reader is directed to Gill, 1982, Ch.ll, for the complete derivation.) The 
governing equation (GE) for LPRW's can be obtained by linearizing the QGPV 
equation. 
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where, 



( ) = d\ )ldy? + d\ )/ayl 

The assumed plane-wave solution for the GE, 

Pq = P(z) exp [ / ( Gt - K ( X cos 6 + y sin 0 ) ) ], 

where the real part, 

Pg = P(z) cos ( Gt - ^x - /y ), 

with 



k = K cos 6, and I = K sin 0, 



is substituted into the GE and simplified, resulting in 



sini Gt - Ax - /y ) { ?(z) [ (A" — 



(j]A/g)] - 



0 

dz ' 



! oP(z) 
dz 



= ). 



which is the GE for the vertical structure. Letting 

-u- = (A- ^ 4- (pA/G), 

gives the resulting form of the equation. 



d 

dz 



( 



1 d? 
dz 



) + H^P = 0 



which is a Sturm-Liouville problem. Applying the rigid (kinematic) boundary 
condition, dpjdz = 0, at z = 0, and at z = -D (where D is the water depth), there will 
be an infinite number of eigenvalues and eigenfunctions, where are the eigenvalues. 
Solving for the dispersion relation (DR) for the nth vertical mode is 



G 



n 









(n = 0,1,2,...). 



The nth mode Rossby radius of deformation can be defined as In 

dimensional form, the barotropic Rossby radius is Xg = where Cg = ( gD)*'^^ 

(wave speed) and /g = 20. sin (tp ) (coriolis parameter, and (p is the latitude). 
Assuming a constant Brunt-Vaisala (N^) profile, the baroclinic Rossby radii can be 
defined as X^^ = ND / nrc/^. By convention, G^ must be greater than or equal to zero; 
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therefore, cos 0 must be less than or equal to zero, (7t/2^0^37i/2), giving a 
westward, with a northward or a southward component, direction of propagation. 



B. THE SYNTHETIC FIELDS 

The synthetic fields were calculated by the following general equation for M 
different waves, each with .V normal modes. Each wave is a suite of progressive modes, 
all having the same and instantaneously in phase at t=0. 



MX. y, z. 
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m= 1 
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where 



F: is the value of the synthetic field at some x, y. z, t, 
is an arbitrarv amniitude function, 

''t>^(z): is a depth function. oosirtz/D), 

is the frequency for the nr/i iiorm.al mode, 
is the total 'vave num.oer. 2“ L. and 
8^: IS the propagation direction. 

The pnysicai constants usea in 'he caicnlations were g = ‘^.Sms'', D = -i km. 
<p = 38. o'", and N = 1.0 x 10‘- s'^ giving /q = 9.07 x 10‘^ s‘% 
P = 1.79 X 10'^^ (ms)‘\ which are representative of a mid-latitude oceanic region as 
in the OPTOMA domains. The limits of the synthetic fields are -32 ^ x, y ^ 288 km, 
to provide flexibility in sampling at the boundaries of the subdomain. The simulated 
sampling subdomain is within the region bounded by 0 ^ x, y ^ 256. 

C. PROCEDURE 

1. Synthetic Field Generation 

Two synthetic fields were generated through use of the general equation 
presented above. Each field was computed at t = 0 ( <7^^ j^t = 0 ) and z = 0 
( Kyz) = 1 ), with = 1 for all waves and all modes. A summary of the specific 

parameters for each field is presented in Table 3. The origin of each wave contributing 
to the field is referenced to the origin of the subdomain, (0, 0). "SYNFLDl" (Figure 
5.1) produced a standing pattern of eddy-like features with a wavelength of ~216 km 
(eddy diameter ~ 108 km). This field was scaled from 10 to 18 ®C, representing a 
nominal SST range as observed in the OPTOMA domains. "SYNFLD2" (Figure 5.2) 
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Figure 5.1 Synthetic field, SYNFLDl and sample subdomain (dashed line), 
representing SST in °C. Contour interval 0.5®C. 
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Figure 5.2 Synthetic field, SYNFLD2 and sample subdomain (dashed line), 
representing MLD in meters. Contour interval 2.5 m. 
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was scaled to represent a mapping of mixed layer depth, with representative values 
from 20 to 70 meters. Again this was based on a similar range of MLD's observed in 
the OPTOMA domains. 



TABLE 3 

SYNTHETIC FIELD SPECIFIC PARAMETERS 
Field: SYNFLDl Number of Wave.v. 2 Number of Modes: 2 



Wave^ 


Origin 


Wavelength (km) 




1 


(-J2, 3 ) 


150.0 


135 


- 


( 11, 106) 


150.0 


225 


id: SYNFLD 


2 'dumber of Waves: 5 


Number of Modes: 


Wave,^ 


vOngin 


Wavelength (km) 


iii'' 


1 


(-23. 52 ) 


200.0 


090 


A 


( J2.-70) 


125.0 


135 




( 6. 23 ) 


1.50.0 


180 




'-65. Z6') 


175.0 




5 


( 74. iiOf 


250.0 


270 



2. Simulated AXBT Surveys 

The simulated survey flights are based on the flight plan employed on 
OPTOMA20P, with flight tracks oriented parallel to the local x axis. 
a. SYNFLDl 

Progressions of six AXBT survey flights, each with a different AXBT spacing were 
conducted in SYNFLDl. Each different spacing was based on the reference spacing of 
Ax = Ay = 28 km, as is used in OPTOMA Program survey flights. Four symmetric 
grids and two "checkerboard' patterns were completed in each progression. 

(1) Progression 1. 

In this progression, four surveys simulated a symmetric x, y 
observation pattern with Ax = Ay = 14, 21, 28, and 42 km. These surveys were 
designated Sl-14, Sl-21, Sl-28, and Sl-42, respectively. The two remaining surveys 
were variations on the 21 and 28 km observations. Every other observation was used 
so that in each row, the samples were offset from each other, as the black squares on a 
checkerboard. These were designated S1-42X and S1-56X. 
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(2) Progression 2. 

In the second progression the same data were used, but a random 
position error, varying between ± 4 km (~2 NM) was introduced to represent aircraft 
navigational error in AXBT locations. These represent one particular realization of the 
random position errors (used in the computations), and were designated 
SM4N,...,S1-56XN. 

b. SYNFLD2 

One survey was completed in this field with Ax = Ay = 28 km. and a A 
km position error. .Again this represents a nonunai OPTOMA survey flight. A 
progression of six objective analysis computations were completed cn the survey data. 
The two versions of the correlation function as presented in Chapter IV were used. 
The input parameters, fraction noise variance and maximum autocorrelation value, 
were variea for each computation. 

3. Application oi Objective Analysis 

For each simulated .AXBT survey in a progression, the objective analysis 
technique was applied as Chapter IV. It is orieily summarized: 

» The autocovariance function and the fitted correlation function are comnuted. 
and confidence imits .at ..-iach oomt on the curve, .cased an the 45% coniidence 
interval for a normal distribution are calculated. 

• From the fitted function, the fraction of the variance due to noise is estimated 
and the zerc*crossing distance is determined. The noise variance, radius of 
influential points, interpolation point spacing, number of influential observations, 
and maximum autocorrelation value are then input into the objective analysis 
program. 

• The objective analysis technique is applied at each interpolation point, and 
quality checks on the value of the autocorrelation matrix determinant and 
maximum autocorrelation values are computed. The objectively analyzed fields 
and their RMS interpolation errors were computed and then mapped. 

4. Comparisons 

For each progression, the correlation (CORR), root mean square error 
(RMS), both the systematic (RMSS) and unsystematic (RMSU) root mean square 
error, and mean absolute error (MAE) were computed between interpolation points of 
each simulated survey and the true field (synthetic field observed at the interpolation 
points). Also each survey was compared with the reference survey 
(Ax = Ay = 28 km). The RMSS error is a measure of the linear bias in the objective 
analysis, while the RMSU error is a measure of the precision of the interpolated field. 
For a complete description of the formulas applied, the reader is directed to Willmott, 
et ai, (1985). 
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D. SUMMARY OF RESULTS 



1. SYNFLDl 

The autocorrelation computations (Table 4) show an increasing estimated 
noise variance with an increase in sample spacing. The worst case was the 42 km 
symmetric surveys, SI -42 and S1-42N, for both progressions 1 and 2. The 
zero-crossing is also much shorter. This could be the result of the larger bin size used 
to compute the autocovariance values as the sample density was more sparse when 
compared to the other surveys. The relative magnitudes of the autocorrelation values 
in each progression are approximately the same. The observation positions, covariance 
function, and fitted correlation function are presented for surveys Si -42, S1-42X, 
S1-56X (Figures 5.3, 5.4, 5.5) and S1-42N, S1-42XN, S1-56XN (Figures 5.6, 5.7, 5.8). 



TABLE J 

AUTOCOVARIANCE AND CORRELATION FUNCTION 
CALCULATIONS, SYNFLDl 



Progression 1 



Survey 


Bin 


Estimated 


Estimated 


Fraction 


Zero 




Size 


Signal 


Noise 


Noise 


Crossing 




(km) 


Variance 


Variance 

"C^ 




Distance 

(km) 


Sl-14 


30 


3.604 


0.316 


0.081 


55.9 


Sl-21 


30 


3.455 


0.483 


0.123 


52.9 


SI -28 


30 


3.081 


0.858 


0.218 


52.0 


SI -42 


45 


1.836 


2.132 


0.537 


39.4 


S1-42X 


30 


3.055 


0.878 


0.223 


59.3 


S1-56X 


40 


2.086 


1.855 


0.471 


56.0 






Progression 2 






Survey 


Bin 


Estimated 


Estimated 


Fraction 


Zero 


Size 


Signal 


Noise 


Noise 


Crossing 




(km) 


Variance 

ocr 


Variance 

o^T 




Distance 

(km) 


S1-14N 


30 


3.551 


0.370 


0.094 


56.2 


S1-21N 


30 


3.584 


0.343 


0.087 


55.1 


S1-28N 


30 


3.262 


0.675 


0.171 


54.2 


S1-42N 


45 


2.115 


1.852 


0.467 


43.6 


S1-42XN 


30 


3.589 


0.343 


0.087 


56.4 


S1-56XN 


40 


2.927 


1.011 


0.256 


60.5 
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Figure 5.3 



O 




0.0 32.0 64.0 96.0 123.0 160.0 192.0 22^.0 256.0 

Distance (km) 




SYN’FLDl, Progression 1. Sample SI -42. (a) Observation positions, 
(b) Autocorrelations (solid) and fitted function (dashed). 

95 % confidence limits (shaded). 
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Figure 5.4 



SYNFLDl, Progression 1, Sample S1-42X. (a) Observation positions, 
(b) Autocorrelations (solid) and fitted function (dashed). 

95 % confidence limits (shaded). 
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Figure 5.5 SYNFLDI, Procression 1, Sample S1-56X. (a) Observation positions, 
(b) Autocorrelations (solid) and fitted function (dashed). 

95 % confidence limits (shaded). 
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Figure 5.6 SYNFLDI, Progression 2 , Sample S1-42N. (a) Observation positions, 
(b) Autocorrelations (solid) and fitted function (dashed). 

95 % confidence limits (shaded). 
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Figure 5.7 SYXFLDI, Progression 2, Sample S1-42XN. (a) Observation positions, 
(b) Autocorrelations (solid) and fitted function (dashed). 

95 % confidence limits (shaded). 



72 




b.O C2.0 64.0 36.0 123.0 160.0 l92.0 224.0 256.0 

Distance ikm) 




Figure 5.8 SYNFLDl, Progression 2, Sample S1-56XN. (a) Observation positions, 
(b) Autocorrelations (solid) and fitted function (dashed). 

95 % confidence limits (shaded). 
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All objective analyses for SYNFLDl were interpolated to a regular 13 by 13 
grid with a spacing of 20 km. This spacing was chosen as it was less than half of the 
zero crossing distance of the reference survey, S1-28N. Input parameters were applied 
differently for each progression (Table 5). An estimated fraction noise variance of 0.05 
and a maximum autocorrelation value of 1.0 were used in Progression 1 to illustrate 
the results of improper specification of the objective analysis input parameters. The 
radius of influence was fixed at 50 km and the zero-crossing distance was 55 km. For 
Progression 2 . the radius of influence for each sur/ey was the computed zero-crossing 
distance and the noise variance was as estimated. The maximum autocorrelation value 
was selected to be 0.7. which limits the distance between selected observations to 
greater than 25 km. Equation 4.2 was the correlation function selected to be 
employed, as in the ADDAS software. 

The objective anaiysis RMS interoolation errors are slightiv larger in 
Progression 2 chan In Progre.ssion !. but the range of values for the determinant is 
decreasea by three orders of magnitude (Table 5). This is a direct resuit of tne 
specification of a smaller maximum autocorrelation value for Progression 2. The 
appearance of the oojectiveiy analyzed flelds 'emains consistent between the two 
progressions tnrough surveys o2-28. Progression 1. and 31-28N. Progression 2. The 
fields become more erratic with the larger observation spacing, SI -42, S1-42X, and 
S1-56X in Progression 1 (Figures 5.9, 5.10, 5.11), while maintaining a consistent 
representation of the field in Progression 2 (Corresponding surveys S1-42N, SI-42XN, 
and Sl-56Xn, Figures 5.12, 5.13, and 5.14, respectively). Again, the cause of the erratic 
nature observed in Progression 1 appears to be the improper specification of the 
maximum autocorrelation value. The fraction noise variance may have a minor effect. 

The comparisons (Table 6) show that the correlation is very high for all 
samples in both progressions. However, samples S1-42N and S1-56XN in Progression 
2 show better correlation than Sl-42 and S1-56X in Progression 1 when "proper" input 
parameters are specified. The magnitudes of the error parameters are generally 
consistent between the two progressions with the exception of the RMS error between 
S1-56X and S1-56XN. The systematic parts are nearly equal whereas the unsystematic 
part is almost five times greater for S1-56X. Since the RMSS is a measure of the 
precision of the interpolated field, and although Progression 2 (S1-56XN) contained 
errors in the positions of the observations, the proper specification, by the user, of the 
input parameters for the objective analysis appears to increase the precision. 



74 



TABLE 5 

OBJECTIVE ANALYSES SUMMARY, SYNFLDl 



Progression 1 



Input Parameters 



Sample 


Radius 


Maximum^ 


Maximum 


Fraction 


of 


Number of 


Autocorrelation 


Noise 




Influence 


Observations 


Value 


Variance 




(km) 






i 


ALL 


50 


'T 

I 


0.7 


0.05 






Computations 


1 

1 


Survey 


RMS 


Vlaximum 


Determinant 


Error 

(%) 


Autocorrelation 


Maximum 


Minimum 


Sl-14 


0.01 


0.918 


0.183E-4 


0.152E-4 


Sl-21 


0.01-0.03 


0.S19 


0.590E-- 


0.294E-4 


SI-28 


0.02-0.15 


0.844 


0.987E-1 


0.I03E-3 


Sl-42 


0.01-0.20 


0.917 


0.850 


0.906E-1 1 


SI-42X 


0.02-0.05 


0.766 


0.650E-1 


0.874E-2 1 


S1-:6X 


0.0 i -0.22 


0.844 


0.937 


0.51^E-2 



Progression 2 



Input Parameters 



Sample Radius 


Maximum 


Maximum 


Fraction 


^ of 


Number of 


Autocorrelation 


Noise 


Influence 


Observations 


Value 


Variance 


(km) 








S1-14N 56.15 


7 


0.7 


0.094 


S1-21N 55.09 


7 


0.7 


0.087 


S1-28N 54.18 


7 


0.7 


0.171 


S1-42N 43.56 


7 


0.7 


0.467 


S1-42NX 56.39 


7 


0.7 


0.087 


S1-56NX 60.45 


7 


0.7 


0.256 




Computations 




Sur. ey RMS 


Maximum 


Determinant 


Error 


Autocorrelation 


Maximum 


Minimum 


(%) 








14N 0.04-0.07 


0.700 


0.107 


0.109E-1 


2 IN 0.04-0.04 


0.696 


0.418 


0.115E-1 


28N 0.07-0.12 


0.699 


0.405 


0.674E-1 


42N 0.31-0.49 


0.446 


4.121 


2.126 


42XN 0.04-0.12 


0.697 


0.174 


0.130E-1 


56XN 0.11-0.26 


0.648 


1.490 


0.505 
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Figure 5.9 SYNFLDl, Progression I, Sample SI -42. 
(a) Objective analvsts. Contour interval 0.5' C. 

(b) RMS interpolation error (%). 
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o 





Figure 5.10 SYNFLDl, Progression 1, Sample S1-42X. 
(a) Objective analysis. Contour interval 0.5°C. 

(b) RMS iriterpoiation error (%). 



77 



o 





Figure 5.11 SYNFLDl, Progression.!, Sample S1-56X. 
(a) Objective analysis. Contour interval Q.f^C. 

(b) R.MS interpolation error (%). 
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o 





Figure 5.12 SYNFLDlj Progression 2, Sample S1-42N. 
(a) Objective analysis. Contour interval 0.5°C. 

(b) RMS interpolation error (%). 
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o 




Distance (km) 



Figure 5.13 SYNFLDl, Progression 2, Sample S1-42XN. 
(a) Objective analysis. Contour interval 0.5°C. 

(b) R.VIS interpolation error (%). 
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Figure 5.14 SYNFLDl, Proeression 2, Sample S1-56XN. 
(a) Objective analysis. Contour interval O.S’C. 

(b) RMS interpolation error (%). 
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TABLE 6 

COMPARISON COMPUTATIONS, SYNFLDI 



Progression 1 



Survey 


CORR 


MAE 


RMS 


RMSS 


RMSU 




(°C) 


CC) 


CO 


CO 1 


SI-14 


0 OOQ 


0.02S 


0.035 


0.022 


0.027 


SI -21 


01998 


0.104 


0.130 


0.040 


0.124 1 


Sl-28 


0.991 


0.197 


0.281 


0.090 


0.267 1 


Si -42 


0.965 


0.360 


0.518 


0.038 


0.517 


S1-42X 


0.997 


0. 165 


0.188 


0.095 


0.162 ! 


S1-56X 


0.379 


0.528 


1.019 


0.379 


0.946 1 

1 






Progression Z 






Survey 


CORR 


MAE 


RMS 


RMSS 


RMSU 






(“C) 


rc) 


CO 


( L) 


SI-I4N 


0.906 


0.14! 


0.180 


0.060 


0.169 1 


Si-21N 


0.998 


f).il3 


0.139 


0.030 


0. 1 36 


S1-2SN 


0.995 


0.166 


0.201 


0.067 


0.189 


S1-42N 


0.994 


0.369 


0.464 


0.041 


0.220 


S1-42XN 


0.996 


0.143 


0.190 


0.086 


0.169 


SI-5dXN 


0.995 


0.321 


0.383 


O.330 


0.194 



2. SYNFLD2 

A single survey, representative of past OPTOMA survey flights, was 
completed. The autocovariance values and the fitted function were computed, resulting 
in a zero crossing distance of 74.9 km and a fraction noise variance of 0.162. The 
objective analysis interpolation points were selected to be at a 32 km spacing in a 
regular grid (Figure 5.15). 

Both versions of the correlation function, as presented in Chapter IV, were 
used with the objective analysis technique for comparison (Figure 5.16). The fitted 
correlation function (Equation 4.1) had a correlation of 0.838 at zero lag as computed 
by the least squares fit (corresponding fraction noise variance of 0.162). In computing 
the fitted function, it was assumed that there was perfect correlation (1.0) at zero lag in 
order to compute an estimate of the fractional noise variance. This value was then 
added to the main diagonal of the autocorrelation matrix. The ADDAS correlation 
function (Equation 4.2) had a correlation of 1.0 at zero lag. 
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Figure 5.15 



SYNFLD2. Sample observations (open symbols) and 
objective analysis interpolation points (dots). 
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Figure 5.16 Correlation functions as used with SYNFLD2. 
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In addition to different correlation functions, the specific input parameters 
that were varied in the objective analysis were the fraction noise variance and the 
maximum autocorrelation value (minimum distance between selected observations). 
These two parameters have been identified as problem areas in the objective analysis 
technique and required further investigation after the surveys on SYNFLDl. Two runs 
of three objective analysis computations were completed. The maximum 
autocorrelation value was changed for the two runs. It was selected from the graphs of 
the respective correlation functions (Figure 5.16), corresponding to distances of 25 icm 
and 35 km. The two correiation functions were applied in the following manner. 

» The fitted function (Ea nation 4.1) was used with the least squares caiculated 
values of the parameters a, d, v. An estimated noise variance, as comoutea bv 
the the least squares fit, of 0.162 (Runs .aI and Bl) was used. 

• The ADDAS function TEquation 4.2) was used in two different wavs; 1) with the 
fraction noise variance selected to be 0.0 (' Runs A2 ana 32), and 2)' with the 
fraction noise variance as computed. 0.i62 tRuns A3 and B3). 

The zero crossing distance remained constant at '4.9 km and the maximum number of 

inriuentiai points was chosen to be seven. Table 7 summarizes the input parameters 

and results of the comparisons. 

All the objectively analyzed adds show a smail range of RMS interpolation 
errors and represent the true fieid veil, with oniy slight differences :n the patterns and 
widely varying error fields. The smallest interpolation errors were observed on runs A2 
and B2 (Figures 5.17 and 5.18, respectively), using the ADDAS function and zero noise 
variance for both the 25 km and 35 km minimum distances. The range of determinants 
for runs A2 and B2 was approximately one order of magnitude smaller than the other 
runs, however. 

In the comparisons, all the runs had a very high correlation with the true field, 
with runs A2 and B2 being the best correlated. The MAE, RMS and RMSU errors for 
runs A2 and B2 were one-half the value of the other runs, also indicating a better 
representation of the true field. Runs A1 and Bl show the highest errors. 

3. Overall 

From the series of calculations presented, it appears that the higher the 
density of the observations, more observations will be available to influence an 
interpolation, and a better objectively analyzed field will result. However, the proper 
specification of the maximum autocorrelation value (minimum distance between 
selected observations) appears to be the key factor in the accuracy of the 
interpolations. Maximum autocorrelation values corresponding to one-third or 
one-half of the zero crossing distance appear sufficient for these data sets and should 



85 



TABLE 7 

CALCULATIONS, SYNFLD2 



Objective Analysis Input Parameters 



Run 


Correlation 


Fraction 


•Vlaximum. 


Corresponding 




Function 


Noise 


Autocorrelation 


Minimum 








Vanance 


Value 


Distance 


A1 


4.1 




0.162 


0.69 


25 km 


A2 






0.0 


0.83 


25 km 


A3 


*n, .2. 




0.162 


0.83 


25 km 


31 


t t 
-f. i 




0.162 


0.56 


35 km 


32 


1 N 
-T. 




0.0 


0.70 


35 km 


33 


-1 




0.162 


0.70 


35 km 






Objective Analysis 


Computations 




Run 




RMS 




Determinant 








Error 




.Maximum Minimum 






i O' 









A1 


5.0-13.0 


0.143 


0.966E-2 


A2 


0.0- 0.4 


0.191E-1 


0.415E-1 


A3 


,5.0- 1-1.0 


0.2569 


0.170E-1 


31 


'.0-16.0 


0.83 / 


0.502E- 1 


32 


0.0- 1.0 


0.^20 


0.112E-2 


33 


7.0-16.0 


0.735 


0.917E-I 



Comparison Computations 



Run 


CORR 


MAE 

(km) 


RMS 

(km) 


RMSS 

(km) 


RMSU 

(km) 


A1 


0.989 


1.168 


1.430 


0.540 


1.324 


A2 


0.997 


0.598 


0.761 


0.111 


0.752 


A3 


0.990 


1.099 


1.358 


0.478 


1.271 


B1 


0.988 


1.299 


1.593 


0.770 


1.394 


B2 


0.996 


0.650 


0.824 


0.179 


0.804 


B3 


0.989 


1.267 


1.526 


0.686 


1.368 



work well in regions of similar variability, using a similar correlation function. The 
noise variance does not appear to be a factor according to the SYNFLD2 results, and 
could also be applied to the SYNFLDl results as well. More work still needs to be 
done with the objective analysis technique and its application to synoptic mapping. 
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Figure 5.17 SYNFLD2, Run A2 (a) Objective analysis (m). 
(b) RMS interpolation error (%). 
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0.0 32.0 64.0 96.0 120.0 160.0 192.0 324.0 356.0 

Distance (km) 




Figure 5.18 SYNFLD2, Run B2 (a) Objective analysis (m). 
(b) RMS interpolation error (%). 
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VI. SUMMARY, RECOMMENDATIONS, CONCLUSIONS 



A. SUMMARY 

With the increasing array of aircraft deployable sensors, systems such as the 
ADDAS, ADAPS, and the Scripps System have proven invaluable in airborne ocean 
sampling and the real-time acquisition and analysis of oceanic data in the research 
communiw. In a matter of hours, aircraft surveys can cover an area that would 
require weeks with a researcn vessel. Synoptic mappings of the oceanic fields obtained 
from AX3T data can provide valuable information for the oceanographer and also 
have a direct application in ASW. 

Proper use of such aata acquisition systems requires a sampling .strategy which 
acquires a representative sample ot'the held that will optimize the information returnea 
with die minimum amount of eifon and time expended. A solid knowledge of the 
dynamics of the oceanic variables m the region must be considered, including both the 
temporal and spatial scales in order to sample the held at a proper interval to 
adequately represent the significant features. ,'\dditionaily, knowledge of the aircraft 
platform and its capabilities, as well as the limitations of the airborne sensors, must 
also be considered. With the data acquired, correct application of the objective 
analysis technique is essential for the accurate representation of the synoptic mappings 
that will be generated. 

B. RECOMMENDATIONS 

1. A Proposed Ship-Aircraft Data Acquisition, Display, and Analysis System 
(SAD AD AS) 

The addition of real-time oceanic data acquisition and processing to a battle 
group can provide a stand-alone capability for the understanding of the oceanic 
environment and the tactical exploitation in ASW. Synoptic mappings in operation 
areas and at strategic transit points would give the warfare commander a "nowcast" of 
the oceanic temperature/sound-speed structure from which appropriate tactics can be 
employed. Furthermore, additional mappings at a later time in the same region can 
provide an assessment of the temporal variability in the region and also of the 
effectiveness of the tactics employed. 
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A system based on the ADDAS, ADAPS, and Scripps System is currently 
possible as these systems are research-proven, both in their hardware and software, and 
would not require a substantial lead time for development for fleet applications. The 
ADDAS is especially feasible as it uses totally off-the-shelf hardware. The ADDAS 
software is based on shipboard data acquisition software developed in the OPTOMA 
Program so a ship/aircraft software package is already available. The software suites 
of both the ADDAS and AD.A.PS are written in Hewlitt-Packard (HP) BASIC and 
could easily be transferred to other HP-BASIC systems. .A software suite could be 
incorporated into the Geophysics Fleet Mission Program Library <'GFMPL), and with 
the addition of a digitizing unit, could provide an interim real-time data acquisition and 
analysis capability until a comprehensive system is available. 
a. System Description 

The main hardware components of the 3ADADAS would be a 
muiti-channei digitizing unit and a multi-tasking computer. The digitizing unit should 
be easily selectable between the array of ur-depioyabie sensors (AXBTs. AXSVs. 
AXSVs. AXCTDs), as well as cxcendaDle sensors such as the X3T and XCTD. With 
a muitichannei system, the simuitaneous processing of two sensors would require a 
computer wath multi-tasking capability. Both ofthe.se components are weil within the 
available technology on the market today. Additionally, a printer and hard-copy 
graphics capability should be included. The HP-9020 and related hardware, currently 
being deployed with the Tactical Environmental Support System (TESS) in the Naval 
Oceanography Community, could provide the hardware capability with only the 
addition of a suitable digitizing unit. However, TESS uses FORTRAN-based software 
so the programs would have to be translated from HP-BASIC. 

The softw'are suite would contain a menu-driven acquisition and storage 
program for each expendable sensor, both air and ship deployable. A program to 
generate the standard Navy and W.MO AXBT- and XBTJJXX messages as the profiles 
are digitized should be included. This would result in less man-hours being used in 
manually selecting the inflection points from a hard-copy profile and provide more 
accurate profiles for archiving at Fleet Numerical Oceanography Center (FNOC). 
Another menu-driven software package would select which field is to be displayed. The 
data would be read from storage and objectively analyzed to produce the synoptic 
mappings. Hard-copy output would then be available to the warfare commander and 
oceanographer for tactical briefing and planning. 
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b. System Employment 

As the name SADADAS implies, the proposed system will be able to be 
deployed in land-based ASW aircraft as well as in ships with deployed air assets. These 
systems are presented below. 

(1) Land-Based Aircraft. 

The aircraft variant of the system would be employed in much the 
same manner as in the OPTOMA Program aerial surveys as presented in Chapter III. 
The hardware and operator would be physically located onboard the aircraft so the 
system must be portaoie and compact. .A system such as this could prove particularly 
attractive in regions of repeated coverage or continuous on-station operations. 
Although there will be some time lag between the initial survey flight and the 
availability of data to subsequent flight crews, synoptic mappings of the oceanic 
variability would be available for tactical planning and could be updated regularly on 
following rlights. A data base could be buiit up and archived at the local .ASW 
Operations Center (ASWOC) for future missions. .A system such as this would 
interface with the TESS rather well. 

(2) Deployed Aircraft Assets. 

A shioboard variant of the .ADDAS could be developed with little 
effort and would provide another facet to real-time data acquisition. The hardware 
and operator would be stationed either in the carrier ASW module or in the Combat 
Information Center (CIC) in destroyers and battleships. AXBTs could be deployed by 
S-3 Viking carrier-based ASW aircraft and the SH-2 Sea Sprite, SH-3 Sea King, and 
SH-60 Seahawk ASW helicopters deployed in frigates and destroyers in the battle 
group. The AXBT audio signal would be transmitted to the ship by the aircraft-ship 
data link. AXBT locations would be provided by either the shipboard air controller or 
by the aircraft onboard navigation and transmitted to the ship. The AXBT audio 
signal would be received from the shipboard sonobuoy lofargram recorder and fed into 
the digitizing unit. Two or more aircraft could be utilized to perform the survey and 
provide larger areal coverage of the region. With the survey complete, near real-time 
synoptic mappings would be available for operational use and tactical decisions. 
Again, this capability could be accomplished with the TESS now being deployed. 

Factors such as the AXBT loadout, data link transmission range, 
navigational accuracy of AXBT positions, and speed and endurance of the helicopters 
need to be investigated more fully to determine the feasiblity of a shipboard data 
acquisition capability and its employment. 
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c. Operationally Applicable Synoptic Fields 

The following fields are presented as tactical decision aids for the warfare 
commander: 

• Sea Surface Temperature. This will indicate the surface manifestations of fronts 
and eddies that may have tactical significance in ASW. 

• Temperature/Sound speed at selected depths. This will determine the depth 
extent of the surface features and provide mappings of tactically significant 
oceanic structure, especially with respect to acoustic propagation. 

• Mixed layer depth. This will provide the warfare commander with the depth of 
the surface duct (if one exists). This has particular application in surface shin 
bow-mounted sonar.„ towed-array sonar, and sonoouov employment, the 
acoustic properties of the mixed laver will' ultimately determine the transmission 
ranges in* the surface duct for bow s'onars. A knowledge of the mixed layer depth 
can provide the towed-arrav operator with an optimum depth to erhploy Ihe 
array. Sonobuoys are limited by the d^th of their hvdrophones. The shallow 
sonobuoys have a hydrophone depth of 90 feet f ~27 meters). With variabilitv in 
the mixed layer, the hydrophone mav be in the layer at one location while below 
the laver at another. 'This will affect the transmission and receiving capabilities 
of the'sonobuoy. 

• Subsurface ducts and sound channels. This will have a similar application as in 
the depth of the surface duct t mixed layer). Subsurface sound channels mav 
affect the the acoustic propagation in the survey region. 

» Vertical sections of temperature/sound speed. Vertical sections through a survev 
region can indicate temperature/sound speed gradients along a specific track from 
the surface to the depth of the survey (normally 300 m). Sound speed profiles 
could also be generated for input to acoustic prediction models* to get the 
transmission loss and acoustic paths along the track. 



C. CONCLUSIONS 

With the capability of research-proven systems and the availablity of off-the-shelf 
hardware and expendable sensors, the development of a SADADAS for the operational 
Navy is currently feasible. The addition of a multichannel digitizing unit and data 
acquisition software to the TESS could provide the Naval Oceanographer and the 
warfare commander with a near-real-time, stand-alone ocean sampling and analysis 
system that can enhance tactical readiness. 

There is an emerging need for ocean sampling and synoptic mapping of the 
mesoscale variability in naval operations, especially in antisubmarine warfare. As 
submarine-radiated noise has been decreasing as advances in quieting technology are 
implemented, the warfare commander will have an ever increasingly difficult task of 
localizing and tracking submarines. Real-time synoptic mappings from AXBT surveys 
can provide the Naval Oceanographer and the warfare commander with a 
three-dimensional representation of the thermal (and hence acoustic) structure of the 
submarine's environment. Standard ASW tactics can be modified and new ones 
developed as the fleet is educated in the variability of the oceans. 
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A most valuable, though unanticipated result of the ADDAS deployment in the 
OPTOMA Program, was that Navy aircrews became aware of the variability of the 
ocean. Usually only one or two AXBTs are deployed in an exercise or mission. These 
AXBTs may not be representative of the ocean structure in the region yet are used to 
compute acoustic parameters which are applied to the sensors. A system such as the 
SADADAS could provide aircrews and shipboard tacticians with a new perspective of 
the ocean, mesoscale variability and its effects on operations. It is time for such a 
system to be developed further and evaluated in actual fleet operations. 
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APPENDIX A 

ADDAS HARDWARE COMPONENTS 

The following components are used in the current version of the ADDAS. Prices 
listed are the aooroximate retail cost as of Julv 1986. Svstem buiid-uo reouires 

■4. 4 . y ^ *.1, 

approximately 40 hours at the GS-7 level. Software testing/ maintenance requires 



about 10 hours at the GS-9 level. 

KP9816 .MC68000 16-bit Microorocessor S5.000.00 

W/HP9121 3.5 inch microfloppy dual disk drive 

Sipoican Mark 9 Front-end Digitizing Unit 5.000.00 

AXBT Circuit Boam 'OO.OO 

HP"^470A Two-pen Plotter 1.200.00 

Avionic Instruments .Vlodei 2A125 -lA Inverter 3,000.00 

VHS Stereo HI-FI Videocassette Recorder SOO.OO 

Tektronix 211 Osciilo.scope 750.00 

Portable Rack-Mount Storage/ Transit Case 500.00 

ARNAV R40 LORAN-C Unit (optional) 400.00 

w/Antenna and cable 

Assorted electronic supplies 200.00 

Computer interface cables 
System wiring 



Total Cost 517,550.00 
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